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{iLi' Abstract 

■ Third QED order hadronic hght-by-hght (LBL) contributions 

aj^^^{M) to the anomalous magnetic moment of the muon a'^"''^ from 
the pole terms of scalar a, ao(980) and pseudoscalar vr'', rj, rj' mesons 
(M) in the framework of the linearized extended Nambu-Jona-Lasinio 
model are evaluated. The off-shell structure of the photon-photon- 
meson vertices is taken into account by means of constituent quark 
triangle loops. The mass of the quark is taken to be niu = rud = nriq = 
(280 =b 20) MeV. The unknown strong coupling constants of 7r^,r],r]' 
and oo mesons with quarks are evaluated in a comparison of the corre- 
sponding theoretical two-photon widths calculated in the framework 
of our approach with experimental ones. The cr-meson coupling con- 
stant is taken to be equal to vrQ-meson coupling constant as it follows 
from the linearized Nambu-Jona-Lasinio model Lagrangian. Then 
one obtains a^^^(7ro)=(8L83±16.50) x lO-^i , a^^^(ry)=(5.62±1.25) x 
10-11, a^^^(r/')=(8.00±L74) x IQ-^^ , a^^^(CT)=(11.67±2.38) x IQ-^^ 
and a^'^'^(ao)=(0.62±0.24) x 10'^^. The total contribution of meson 
poles in LBL is a^^^(M)=(107.74 ± 16.81) x lO^^^. 
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1 Introduction 



The muon is described by the Dirac equation and its magnetic moment is 
related to the spin by means of the expression 

f^ = 9(7^]s (1) 

where the value of the gyromagnetic ratio g is predicted (in the absence of 
the Pauli term) to be exactly 2. 

In fact, however, the interactions existing in nature modify g to be ex- 
ceeding the value 2 because of the emission and absorption of virtual pho- 
tons (electromagnetic effects), intermediate vector and Higgs bosons (weak 
interaction effects) and the vacuum polarization into virtual hadronic states 
(strong interaction effects). 

In order to describe this modification of g theoretically, the magnetic 
anomaly was introduced by the relation 

= ^= (2) 



where a = 1/137.03599976(50) fTj] is the fine structure constant. 

Here we would like to make a note that from all three charged leptons 
(e~, r^) the muon magnetic anomaly is the most interesting object for 
theoretical investigations due to the following reasons: 

i) it is one of the best measured quantities in physics P 

a^^P = (116592020 ± 160) x 10"^^ (3) 

a) its accurate theoretical evaluation provides an extremely clean test of 
"Electroweak theory" and may give hints on possible deviations from 
Standard Model (SM) 

in) the new measurement |^ in BNL is expected to be performed with a 
definitive accuracy 

Aa'^^^'P = ±40 X 10"" (4) 

i.e. it is aimed at obtaining a factor 4 in a precision above that of the 
last E-821 measurements (H). 
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At the aimed level of the precision (^) a sensibility will already exist to 
contributions 0, |^ 

^(2,3)«,eafc _ ^^53 ± 4) X 10"", (5) 

arising from single- and two-loop weak interaction diagrams. 

However, as we have mentioned above, the muon magnetic anomaly may 
also contain contributions from a new physics, which can be revealed, we 
hope, in a comparison of a^^'^ with an accurate theoretical evaluation of a^l^. 
While QED and weak interaction contributions to a^^ seem to be estimated 
very reliably, there is still opened a door for improvements in hadronic con- 
tributions. 

As the most critical from all hadronic contributions are the light-by-light 
(LBL) meson pole terms [Q, we recalculate the third-order hadronic LBL 
contributions to the anomalous magnetic moment of the muon al^'^ from the 
pole terms of the scalar a, ao and pseudoscalar 7r°, 77, ry' mesons (M) in 
the framework of the linearized extended Nambu-Jona-Lasinio model. The 
reason for the latter are predictions of series of recent papers 



«f^(vro) 



-52 ± 18) X 10"^^ [7] 

-92 ±32) X 10"^^ [8] 

-79.2 ± 15.4) X 10"^^ [9] (6) 

+83 ± 12) X 10"" [10] 

+58 ± 10) X 10"" [10,11] 



which differ not only in the magnitude, but even in the sign. Moreover, in 
these papers only the pseudoscalar pole contributions were considered. In 
our paper we include the scalar meson (cr, oq) pole contributions as well. 

The crucial point is description of transition form factor 7* — > M7*. 
Current methods used for this aim are the ChPT and the vector-meson- 
dominance (VMD) model. 

In this paper the corresponding transition form factors by the constituent 
quark triangle loops with colourless and flavourless quarks with charge equal 
to the electron one are represented. (An application of a similar modified 
constituent quark triangle loop model for a prediction of the pion electro- 
magnetic form factor behavior can be found in [|12|, where also a comparison 



with the naive VMD model prediction is carried out.) The mass of the quark 
in the triangle loop is always taken to be m„ = nid = niq = (280 ± 20) MeV 
||13|| , which was determined in the framework of the chiral quark model of 
the Nambu-Jona-Lasinio type by exploiting the experimental values of the 
pion decay constant, the p-meson decay into two-pions constant, the masses 
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of pion and kaon and the mass difference of rj and 1]' mesons. The unknown 
strong couphng constants of 7r°, rj, rj' and mesons with quarks are evalu- 
ated in a comparison of the corresponding theoretical two-photon widths with 
experimental ones. The a-meson coupling constant is taken to be equal to 
7r*^-meson coupling constant as it follows from the corresponding Lagrangian. 
The (T-meson mass is taken to be mo-=(496 ± 47) MeV as an average of the 



values recently obtained experimentally from the decay 7r^7r+7r+ [14 

and excited T decay processes. 

As a result we present explicit formulas for a^^^{M) (M = vr°, r], rj', a, oq) 
in terms of Feynman parametric integrals of 10-dimensional order, which 
subsequently are calculated by MIKOR method [T^. 

The paper is organized as follows. In the next section all definitions and 
a derivation of basic relations is presented. More detail can be found in 
Appendix A and in Appendix B. The last section to numerical results and 
discussion is devoted. 



2 Meson pole terms in light— by— light contri- 
butions to a^^^ 

The third QED order hadronic LBL scattering contributions to a^"''^ are gen- 
erally represented by the diagram in Fig. |l|^ , which contains a class of 
pseudoscalar mesons (vr''^, K"^ and also Kgj^) square loop diagrams (Fig. ^B), 
a class of quark square loop diagrams (Fig. |l|C) and scalar a, ao(980) and 
pseudoscalar tt^, r], t]' meson pole diagrams (Fig. where the off-shell struc- 
ture of the photon-photon-meson vertices is taken into account by means of 
flavourless and colourless constituent quark triangle loops. Here the interac- 
tion of mesons (M) with quarks is described by the linearized Nambu-Jona- 
Lasinio type Lagrangian 

^qqM = guqix) [uix) + ivr(x)75] g(x), (7) 

with unknown strong coupling constant Qm- The latter for vr*^, r], rj' and 
oo mesons are evaluated by a comparison of the corresponding theoretical 
two-photon widths with the experimentally determined values [|I|. 

For the on-mass shell scalar {S) and pseudoscalar (P) mesons decays one 
can write the following matrix elements 

M{S{p) ^ 7(fci) + 7(^2)) = --^K{ms/m,){k,k2g,u - (8) 

-kivk2f,)ei^{ki)e2^,{k2), 
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ABC 

Figure 1: Third order hadronic light-by-light scattering contribution to a|j'^'^ 
(A) and class of pseudoscalar meson square loop diagrams (B) and quark 
square loop diagrams (C) contributing to (A). 



M{Pip) ^ j{h) + 7(^2)) = ^Amp/mg){hhe,e2), (9) 

with 



1 ^-x _ . _ 1 







1 — xyz 



and 



= m^, kf = kl = 0, (abed) = eap-yscfh^ . 

They lead to the theoretical two-photon widths (for more detail see Ap- 
pendix A) 

r? = (11) 



r? = ^^^^^(".p/m,). (12) 



647r^m 



dependent on the scalar gs and pseudoscalar gp meson coupling constants 
with quarks, respectively. 

Taking the meson masses with errors and the following world averaged 
values of the two-photon widths from the last Review of Particle Physics 
r;;:^: = (O.OOS ± O) keV, r;^^ = (0.464 ± 0.044) keV, VJ^ = (4.282 ± 0.339) keV 
and r2;j = (0.24 ± 0.08) keV, the quark mass value rrig = (280 ± 20) MeV, 
one finds the meson coupling constants with quarks as follows 

glo = 9.120 ± 1.305, g^ = 2.220 ± 0.378, (13) 
g^ = 6.708 ± 1.096, gl^ = 0.757 ± 0.277. 
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Since there is no experimental result on the cr-meson decay into two- 
photons up to now, we identify with g'^o as it follows from the Lagrangian 
(7). 

Now we are ready to evaluate LBL meson pole terms contributions to 



„had 



Usually the fermion-photon vertex function is written in the form 

T, = u{p2)VM)u{p^), (14) 
= 7M^i(g') + 4^[g,7M]^2(?'), ? = P2 - Pi, 

with on mass-shell muons, i.e. p\ — p^ — {fn is the muon mass). Then 
the expression for the muon anomalous magnetic moment is 

a, = ^^^ = ^2(0) (15) 

3 1 

lim —Sp{p2 + m)V^{q){pi + m)V^{q), 



with a projection operator 



The gauge invariant set of Feynman diagrams containing the LBL scat- 
tering block with meson pole intermediate states is drawn in Fig. 2. 



Their contributions to a}}°^ are 



(16) 



lim ^Sp{p2 + + m)V^{q), 



with 



yS+P 



/ 



2^2 



2 \a 



k'l — nip kf —m 



klkl{k^ - k2y 
O = 20i + O2, 



, (17) 



QvXa ^ P2-k2 + m pi + fci - fca + m 
^ (P2 - ^2)^ - rn'^ (Pi + ki — k2Y — w? 



where 



T. 



(P) ^4a^ 



{ljiykiq){\ak2ki)fi{kl)f2{ki, /C2), 



(18) 



with 



/i(^i) 



1 l-xi 



^ — a;(l — a;)/c' 



dxi J -^dx2, (19) 



d = mq- klx3{l - xs) - klx2{l - X2) + 2/^1/^2X2X3, 

X3 = 1 — Xi — X2. 

Here we would like to note that the functions fi{kf), f2{ki,k2) at the limit 
of an infinite value of the quark mass, niq — > 00, tend to one. 
The expression for T^^^ has a similar form 



HfXa 



4a^ 



{Qnukiq - ki^q^)fi{kl)fxa{h, ^^2), 

m^(l — 2x(l — x))dx 



(20) 



— x(l — x)A;i 

1 l—xi 2 



2 ' 



/r in 
dxi J {{I - 4:XiX2) [gxak2{ki - k2)- 



k2x{h - k2)a] + 2(/i;i - A;2)^xi(2xi - l)gxa} dx2. 



Usual procedure of performing the loop-momenta integration using the 
Feynman parameters method (see Appendix B) leads to 



where 



and 



a^^^(iV/) = -RMFM{mM,mg) 



R 



M 



647r^m^ 



Fp(mp, rriq) = 2M„ + M^, Fs{ms, rUg) = 2Na + N^. 
The Mj, Ni are given by the expressions 



M. 

N, 



' A^ 

' A^ 

, r5 





1 , 1 

ax 





xi 1 Zl 



/ (iFf, = I — I dxi I ^ / / ^2:2 / 

y J XX J J X3X3 J J J 



"^"^ 



dzi J dz2 / a2Ci^Ci^d app dp, 




a; = 1 — x, 2;i = 1 — ^1, r 



(21) 
(22) 



(23) 



with dra,b the 10-dimensional integrals on Feynman parameters 

1^1 ^ 1 Zl 1 

/ (iFa = t'^ — dxi / — ^ / dzi I dz2 / a2a\d^adypp'^dp, (24) 

t/ J XX J J X^X^ J J J 





, d^a = daida2da3, d'^a = daida2da3da4. 



m 



The explicit expressions for polynomials A, B are 



8yp 



2a\ + 2riai{ai + hi) -rlaihi - Srf, 



4 + 12ri + 4ai6i^, 

z 

^[4 - rla2h2{h2 - 1)] - yr^as + ^rl{l - 262), 
4a2 - ^ + r2(862 - 4) + — 0262(^2 - 1), 



with 



^6 



+ 



+ 



yp 

z 



-8al - 16riai(ai6i - ^ai + ^61) 



8rl{aibi ai H — bi )bi 



+ — 



yp 

z 



3 

32 



di{xi,X2) 



(16/3)a^(ai + 1) + —n{bi - l)(ai + l)ai 



+ rffia, + l)ib,-ir 



3 

0^2 (Xi), 



16. 



(8ai + y)ci2(2;i) - (4 + 12ai)cii(a;i, X2) 



8yp 



yp, 



ia, + l){b,-lfd2{xi) 



+ —{Uaibi - 8ai + 861 - 8)bidi{xi, X2), 
z 

P f [16 2/ -,\ 32 , ^, 
= -\ ^a2(«2 - 1) + ^r2a2(a2 - 1)(02 - 1) 
(7 I L o o 

+ yr2(a2 - 1)(62 - 1)' 



16 a 16 o, ^, 

-—al + —rlb2{b2 - 1) 



--^62(62-1) 

(T 

8a2-^-^(a2-l)(62-ir 



^2(2^1) 

c?i(a;i,a;2)| , 

di{xi,X2) 

d2{xi) 



Oi = q;3Q;2, bi — q;2Q;37i, 02 = a4{a3a2 + 0:3), 
62 = tt4(tt37i + "3), n = piy - Civ), r2 = -C2P, 
di{xi, X2) ^1- 4x1X2, d2{xi) = 2x1(2x1 - 1). 



Quantities Ci^2, 7i, z, Di, D2 and a are given in Appendix B. 
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3 Results and discussion 

The result for the meson pole contributions of LBL type has the form 

a^^'(M) = ^(7r°) + ^(ry) + af\v') + a^^^\a) + ^(ao) (27) 



where 



a 



= -Rr^'Fr^irrir^^^mg) (28) 



«^^^(f^) = -RaFJm^,m„ 



o-i^^M = -RaoFaoimao^rriq). 

Now, taking the values (13) of g\,i determined in Section 2, the mass 
of the quark uiq = (280 ± 20) MeV ||l3l, the mass of the cr-meson mg- = 
(496 ± 47) MeV |T5| , the mass of the muon and the corresponding scalar 
and pseudoscalar mesons from the Review of Particle Physics one finds 

aj^-^(7r°) = (81.83 ± 16.50) X 10"" 

af^^^{ri) = (5.62 ± 1.25) X 10"" 

a^^^{r]') = (8.00 ± 1.74) X 10"" (29) 

a^^^(a) = (11.67 ± 2.38) X 10"" 

aj^^(ao) = (0.62 ± 0.24) X 10"". 

Here we would like to note that our value for the vr" is in the framework of 



error bars consistent with 10, 11 



So, the total contribution of meson poles in LBL is 

aJ^-^(M) = (107.74 ± 16.81) x 10~", (30) 

where the resultant error is the addition in quadrature of all partial errors of 
(30). 

Together with the contributions of the pseudoscalar meson (vr^, K^, Kg ^) 
square loops and constituent quark square loops |0, || it gives 

al^^itotal) = (111.20 ± 16.81) x 10"". (31) 

The others 3-loop hadronic contributions derived from the hadronic vac- 



uum polarizations (VP) were most recently evaluated by Krause JTT 



a(3)^^ = (-101±6) X 10"". (32) 
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Then the total 3-loop hadronic correction from (32) and (33) is 

^{3)w ^ a'f'itotal) + a(f)^^ = (10.20 ± 17.28) x 10"" (33) 

where the errors have been again added in quadratures. 

If we take into account the most recent evaluation [113] of the lowest- 
order hadronic vacuum-polarization contribution to the anomalous magnetic 
moment of the muon 

^{2)had ^ ^-^Q21 ± 76) X 10"" (34) 

(it is interesting to notice that this result almost coincides in the central 
value and in the error as well with the result [|l9l of one of the authors (S.D.) 
of this paper obtained 12 years ago) the pure QED contribution up to 8th 
order [pO[ 

^QED ^ (116584705.7 ± 2.9) x 10"" (35) 

and the single- and two-loop weak interaction contribution finally one 
gets the SM theoretical prediction of the muon anomalous magnetic moment 
value to be 

af = (116591888.9 ± 78.1) x 10"". (36) 
Comparing this theoretical result with experimental one (3) one finds 



a 



which implies very good consistency of the SM prediction for the anomalous 
magnetic moment of the muon with experiment. 



Appendix A 

We put here the details of calculations of the two-photon widths in the 
framework of the linearized Nambu-Jona-Lasinio model. 

The decay matrix element of a pseudoscalar meson P into two photons 
has a form (we take into account two directions of fermion line in the triangle 
fermion loop) 

M(P(p) ^ = / (A.1) 

with 

5*1 = -^Tr{k -k2 + ■m„)e2{k + m„)ei(A; + fci + m„)75 = i{k2kieie2) (A.2) 
ai = {k + kiY - rrig, 02 = {ki - ^2)^ - m^, = k'^ - m^, 
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where we use the definitions 75 = i7°7V^7^, \Spi^nul\la = i^ixvXa, £0123 = 
+1. 

Joining the denominators and performing the loop momentum integration 
one obtains 

^ 2gpa r ci^X(5(E - 1)^^(^2^16162) 
^ Trniq J — XiX2P^ — XiXa/c^ — X2X^k'l 

This expression is used for the construction of f2{ki, k-i) in (19). 
For the case of real particles we have 

Ml'' = ^{k^k2e^e2)J{mp/m,) (A.4) 

with J{z) given by (10). 

Note that for the case of T^'-meson the matrix element becomes complex. 
Numerically one obtains 

J^(m^o/mg) = 1.04, J2(m^/mg) = 3.73, J(m^//mg) ^ = 2.12 . (A.5) 

For the matrix element of a decay of the scalar meson S into two photons 
one can write 

M{S{p) ^ 7(A;i, ei)7(A;2, 62)) = e^^e^.M^^^ (A.6) 

where 

TT J 010203 
Ti''' = j^SpCk -k2 + mg)^^{k + mq)^^{k + ^2 + rUq). 

This quantity suffers from the ultraviolet divergences. It is the gauge in- 
variance requirement which provides the regularization. General form of 
M'^^{ki,k2), which satisfies the current conservation condition 

M'^^(A;i,A;2)[A;i^;A;2.] = 0, (A.8) 

has the form 

M'"'{ki,k2) = A{g'"'kik2- k'^k!^) (A.9) 
+B{klk^ - kik2k\){klk^ - kik2k!^). 

Usual procedure of joining the denominators, the calculation of the trace, a 
performation of the loop momentum integration and an application of the 
gauge conditions to the muon block, leads to a final result. 
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For the case of real particles the structure B disappears and we have 

^(5^77) = e'l{h)e'^{k2){g^.hk2 - klk!^)^^K{ms/m,), (A.IO) 



nniq 



which leads to the expression (11) for the width with K{z) given by (10). 
Numerically 

\K{mao/mq)f = 0.96. 

Appendix B 

For a performance of the loop momenta integration we use Feynman trick of 
"joining" the denominators 



1 (n + m-1)! } x'^-^il-xY'-^dx 



^mij,n (n - l)!(m - 1)! y (ax + 6(1 - a;))'"+" 
d^k (A;2)"^ (^+i)!(^_^_3)! 



(B.l) 



and the relation 

J in^ (F -d)-~^ ^' {n - l)\d"-2 ,ri>m + Z. [t5.Z) 



For this aim we present the denominators containing k2 in the form 
1 1 



{5}' {6}, 



B.3) 



-a;3a;3{l}{2}{3}{4} 
where 

{1} = kl {2} = kl - 2p2k2, {3} = {k2 - k,)\ {4} = rf/(-X3X3) 

{5} = kl - 2piki; {6} = {pi + ki- k2f - x = l-x. 

Using Feynman parameter a\ to join the denominators {3}, {4} and the 
parameter a2 for joining the result with {2} we have 

a2[ai{4} + q;i{3}] + a2{2} = kl - 2k2a + a2Ai, a = 0:271^1 + 02^2, (B.4) 

X2X2kl _ _ X2 _ ^ 

Z\i = a20ii h a20i2Ki-i ^\ — ai -\- a\ — ^ a — i. — a. 

X^X^ X3 

For the diagram in Fig. 2 A we use the parameter 0:3 for the joining of the 
denominator {1} with the above result 

{k2-hf-zd^, (B.5) 

, , 2 /- , OiiX2X2. 

X3X3 

di^ kl + 2cikipi + A2, 

Ot2(^'iot2''y\ . _ 1 r 2 2-2 I 2'^l'^2tt3 



^z^3^z / 1 A 
Cl = , ^2 — 

z z 
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m a^a2 + m 



X3X3 



For the diagram in Fig. 2B we use the parameter for the joining of the 
denominator (B.3) with {6} and parameter to join {1} with the result 

{k2 - 62)' - (Jd2 (B.6) 



where 



a = al^as'ji + 0:3)^ — a4^as — a2CfiOii ^ai + a 

62 = 04 [(as + ai7i)/ci + (^3 + Q;3«2)Pi] , 
d2 = kl + 2c2Piki + 5, 



C2 = -a4 [q;4(q;37i + Ois){a3a2 + as) - as 
a 

1 



X2X2\ 
^XsXsJ ' 



5 = - 
a L 



m a^{asa2 + as) + 



aia2asa4 2 



XsXs 



To perform the ki integration we use the parameters zi,Z2 to join the 
denominators containing kf 



1 



1 



XX kl{kl - M^){kl - Ml) XX 



dzi 



with 



dZ2 



(B.7) 



XX 



ds^kf- ziMt - Z2M\ M( 

For the diagram in Fig. 2 B we use parameter p to join the denominators 
d2,d^ 

pd2 + pds = (/ci - q2f - m^D2] (B.8) 
?2 = -PC2P1, m^D2 = w?{c2pY - 5p + p{ziMl + Z2M'^). 
For the diagram in Fig. 2 ^ we use parameter y to join di with {5} 

ydi + y(5) = - 2piki{y - ay) + yAi (B.9) 

and parameter p to join the result with d^. As a result one gets 

{ki - qif - Dirri^; qi = p{y - Ciy)pi, (B.IO) 
m^Di = m^p^{y - yc^)^ + p{ziMl + Z2M'^) - pyAi. 
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